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ABSTRACT 


an  on-off  control  provides  a  minimum  time  response  for 
missile  control.  For  application  in  missile  control  systems, 
it  is  wasteful  of  control  effort  (due  to  chatter)  to  use  a 
ideal  relay.  Hence   it      is   necessary      to    modify      the    ideal 

relay  into      a   saturating      linear   control.  The   result      was 

almost  the   same    to    that   of    using    the    ideal   relay. 
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I.    INTRODUCTION 

Guided  missiles  are  classified  into  four  broad  catego- 
ries, depending  on  launch  and  target  position  characteris- 
tics. These  categories  are  (1)  .  air  -to  -  air  (2).  air  - 
to  -  ground  (3)  .  surface  -  to  -  air  and  (4)  .  surface  -to  - 
surface.  Each  category  of  the  abov=  will  employ  one  or  more 
of  the  following  guidance  schemes;  programmed  command, 
liae-of-sight,  lead-angla,  proportional  navigation  homing 
and  inertial.  The  beam  rider  guidance  is  included  in  the 
line-of-sight   guidance.  A   number   of      missiles    also      use  a 

combination  of  these  methods.  For  axample,  the  initial  part 
of  the  missila  trajectory  may  use  programmed  guidance  while 
the    terminal  phase    may  use    beam-ridsr. 

This  thesis  discusses  the  surf ace-tc-air  missile 
controlled  by  on-off,  thrust  vector,  control.  Consideration 
was  given  to  determine  the  effects  of  the  two-level  relay 
and  the   saturation    linear  control.  In   order  to    verify   the 

results,  it  was  tested  by  using  the  type  of  control  for 
three   different    types  of   missile-target   scenarios: 

(1).      LOS   command  against    non-maneuvering   target 

(2) .      LOS   command  against    maneuvering   target 

(3)  .      Pseudo-LOS    command  against    non-maneuvering   target. 

In   chapter    2,      a    discussion      of   a   line-of-sight    guidance 

was    presented      and    a      practical  example    of      it   was      shown   in 

chapter      3.        The     general      concept      of   on-off     control      was 

described   in  chapter   4.      The   simulation    results   of   the    basic 

command   to   line-of-sight   against      both   a   non-maneuvering  and 

a  maneuvering     target  were    shown   in      chapters   5   and      that   of 

pseudo-LOS   case    was    in  chapter   6.      Finally,      a   discussion   of 

two-level  relay    and    saturating      linear   control   was    presented 

in    chapter    7.         A   table   of    variables    which    were    used   in  this 


thesis   is   shown    in    the     Appendix   A.      Computer   simulation   was 
accomplished  using    Digital    Simulation   Language,    DSL. 
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II.    QVEBVIEW    OF   LINE-OF-SIGBT    GUIDANCE    CONTROL 

A  LOS  system  can  be  called  a  "3-point"  guidance  system 
since  there  is  one  point  which  defines  the  tracker,  another 
the  target  and  a  third  which  defines  the  position  of  the 
missile.  The  object  of  the  guidance  system  is  to  constrain 
the  missile  to  lie  as  nearly  as  possible  on  the  line  joining 
the  tracker  and  the  target  called  the  Line  Of  Sight  (LOS). 
The  concept  is  simple  and  can  be  implemented  in  many  ways; 
perhaps  it  is  this  apparent  simplicity  which  explains  why 
many  of  the  guided  weapon  systems  as  yet  designed  are  LOS 
system. 

Consider  a  target  flying  straight  and  at  constant  speed, 
and  a  missile  flying  at  a  different  angle  but  constant 
speed,  having  been  launched  when  the  target  occupies  a  posi- 
tion   TO     (see  Figure    2.1). 

After  intervals  of  time  of  1,2,3  etc  seconds  the  LOS  is 
shown  as  OT1,OT2,OT3  etc.  Since  the  missile  ideally  always 
lies  on  these  lines  the  flight  path  will  be  a  curved  one, 
for  an  approaching  target,  the  curvature  becomes  increas- 
ingly severe  towards  the  end  of  the  engagement.  We  note  that 
the  tangent  to  the  flight  path  at  any  one  point  defines  the 
instantaneous  direction  of  the  missile  velocity.  It  is  seen 
that  the  missile  velocity  vector  will,  in  general,  not  be 
directed  along  the  LOS;  towards  the  end  of  the  engagement  it 
may    be   at  a   considerable  angle    to   it    [Bef.    1], 

In  an  actual  situation  the  guidance  signals  transmitted 
to  the  missile  are  the  demanded  lateral  accelerations 
(LATAX)  in  two  axes  at  the  right  angles  to  the  beam.  These 
demands  are  resolved  into  missile  axes  within  the  missile. 
An  error  compensation  term  endeavouring  to  keep  the  error 
off    the   beam    (  &mj.)    equal  to    zero. 
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Figure   2.1        Hissile  Target   Encounter    with   LOS   Guidance 

A  basic  geometry  and  a  simplified  guidance  loop  are 
shown   in   Figure    2.2    and   Figure   2.3. 

Suppose  that  the  cross  range  error  (CRE)  of  Figure  2.2 
can  be  measured  either  directly  or  by  means  of  the  angular 
difference  between  OT  and  ON,  together  with  some  knowledge 
of   missile   range    (Rm) ,   then 


CRE     =    RMtft  -tfm     ) 


(2.1) 


If  this  error  off  the  beam  is  used  is  an  acceleration  demand 
U,  it  needs  some  damping  so  that  good  response  characteris- 
tics   are  obtained.    A    dynamic   equation   of    the    form 


CRE    ■    G1    (CRE)     ♦      G2     (CRE) 


(2.2) 


needs  to  be  satisfied,  where  G1  and  G2  are  constants.  This 
necessity  leads  immediately  to  the  consideration  of  a 
filtered   error.      In    the   presence    of   noise   on   the    sight-line. 
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and  hence  on  the  cress  range  error,  CRE,  such  a  filtar 
design  is  not  simple  and  becomes  a  compromise  between 
requirements  for  smoothing  the  noise  and  giving  an  adequate 
response  to  a  demand.  Modern  techniques  allow  filters  to  be 
designed  statistically  if  some  knowledge  of  the  noise  char- 
acteristics is  available  or  can  be  assumed.  Figure  2.3  shows 
the  position  of  such  a  filter  F(s)  in  the  guidance  loop.  It 
includes   a   gain    G,    and  the    acceleration   demand  is 

0      =   F(s)     Rm    (    <5t  "  6m    )  <2-3> 

The  missile  transfer  function  is  represented  by  A(s)  and 
when  the  achieved  acceleration  is  doubly  integrated  and 
divided  by  Rm  it  represents  a  new  measure  of  the  missile 
beam  angle  (tfm)  ,  thus  closing  the  loop  when  differenced  with 
the    target    beam    angle    (dt)  . 

While  this  concept  is  simply  a  LOS  or  beam  riding 
guidance  situation  it  is  by  no  means  as  clear  in  homing  how 
a  guidance  law  can  be  devised  in  the  absence  of  information 
on   missile    and   target   positions  [Ref.    2]. 
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Ill-     TYPICAL    ENGAGEMENT    SEQUENCE 

In  order  to  provide  a  "vehicle"  through  which  to  better 
understand  the  basic  aspects  of  command  to  line-of-sight 
guidance  methodology,  the  engagemaat  seguence  of  a  short- 
range,  air-defense,  missile  system  is  iescribed.  The  Roland 
system  was  selected  because  the  general  operational  aspects 
of  the  system  are  available  at  the  unclassified  level 
[Ref.   *]. 

The  entry  of  one  or  more  aerial  targets  into  the  range 
of  the  search  radar  is  indicated  to  the  Roland  vehicle 
commander  by  an  audible  tone.  At  the  same  time,  a  synthetic 
display  of  the  targets  appears  on  a  screen  to  give  the 
commmander  all  the  information  needed  to  select  the  most 
threatening  target.  The  screen  images  are  different  for 
friendly  and  enemy  targets.  Also,  the  entry  of  the  target 
into  the  missile  envelope,  utilizing  target  advanced-range 
computations,  is  indicated  by  a  change  in  the  display.  With 
the  search  antenna  raised  and  the  search  radar  activated, 
target  acquisition  is  possible  even  when  the  vehicle  is  in 
motion. 

There  are  three  modes  of  identification,  friend  or  foe 
(IFF)       interrogation:  automatic,      manual,        and   automatic 

within  a   given    range. 

When  the  commander  has  recognized  a  target  as  hostile 
and  decided  to  engage  it,  he  places  a  cursor  over  the  screen 
image,.  This  automatically  brings  the  turret  to  bear  and 
tracking  can  commence  in  either  the  "radar"  or  "optical" 
modes. 

In  the  "radar"  mode,  the  tracking  radar  automatically 
accepts  target  designation  from  the  search  radar,  searches 
for,    locks   onto,    and   tracks    the   target. 
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In  the  "optical"  mode,  the  aimer  searches  for  the  target 
in  elevation  with  an  optical  sight.  To  aid  him  an  electronic 
instrument  displays  the  maximum  theoretical  elevation  for 
the  search.  ahen  the  aimer  has  acquired  the  target  in  his 
cross-hairs,  he  keeps  the  target  in  his  sight  by  manipu- 
lating a  control  stick.  This  control  keeps  the  target  prop- 
erly positioned  by  moving  the  turret  in  azimuth  and 
swivelling    a  mirror    in  elevation. 

As  soon  as  the  commander  confirms  that  the  target  is 
within  missile  range,  he  initiates  the  firing  seguence  in 
the  "radar"  mode,  or  authorizes  "optical"  mode  firing 
through  a  command  displayed  in  the  aimer's  sight.  The  aimer, 
then,    can   initiate   the  firing    sequence. 

The  missile  is  guided  by  a  command  to  line-of-sight 
technique.  This  means  that  the  target  is  tracked  optically 
or  by  radar  and  the  deviation  of  the  missile  from  this  line 
of  sight  is  determined  and  corrected  by  a  guidance  command. 
The  commander  may  switch  from  "radar"  to  "optical"  and  back 
again,    as   desired,    even   after   the   missile    has   been    launched. 

Target  tracking  and  determination  of  the  missile's  devi- 
ation from  the  line  of  sight  are  different  for  each  mode.  In 
the  "radar"  mode,  the  guidance  radar  has  two  receiving 
channels.  One  is  used  for  target  tracking  and  the  other  is 
used  to  locate  the  missile  in  the  radar  lobe  through  recep- 
tion of  the  missile's  radio  frequency  beacons.  By  comparing 
these  angles,  an  error  between  the  missile  and  the  target 
line  of  sight  can  be  determined.  In  the  "optical"  mode,  a 
biaxially-stabilized  mirror  is  manually  controlled  to  keep 
the  target  vertically  in  the  aimer's  sight  and  the  turret  is 
rotated  to  the  azimuth  of  the  target  line  of  sight.  An 
infrared  goniometer  is  mounted  to  provide  misile  angle  from 
the  tracker  by  following  flares  mounted  on  the  rear  of  the 
missile.  Then,  a  deviation  of  the  missile  angle  from  the 
target   line  of   sight    can   be    determined. 
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Two  groups  of  signals  are  introduced  into  the  command 
computer:  the  velocity  of  the  line  of  sight  in  azimuth  and 
elevation,  and  the  deviation  of  the  missile  from  xhe  line  of 
sight  in  azimuth  and  elevation.  Based  upon  data  from  the 
line-of- sight  movement  and  the  angular  deviations  of  the 
missile,   the  necessary  guidance   signals    are  calculated. 

The  guidance  signals  are  relayed  to  the  missile  by  a 
transmitter  with  highly  directional  characteristics.  The 
command-transmitting  antenna  is  slaved  to  the  missile  angle 
in  both  azimuth  and  elevation.  It,  therefore,  is  trained  on 
the    missile   continuously. 

The  side  forces  required  for  missile  course  corrections 
are  produced  through  deflection  of  the  exhaust  jet  of  the 
sustainer  motor  by  spoilers  at  the  rear  of  the  missile 
(thrust-vector    control)  . 

When  the  missile  reaches  the  point  of  impact  with  the 
target,  the  warhead  is  detonated  by  either  percussion, 
contact  fuse  or  the  radio- frequency ,  proximity  fuse.  The 
warhead  consists  of 'a  radial-effect,  multiple-fragmentation 
charge. 

Figure  3.1  presents  an  operational  schematic  of  the 
basic   Roland  missile   system    operation. 

The  computer  simulations  contained  herein  are  generic  in 
nature  within  the  command  to  line-of-sight  guided-raissile 
type  and  have  only  reasonable  estimates  of  missile  capabil- 
ities introduced.  This  ensures  unclassified  results.  At  the 
same  time,  the  simulations  are  of  sufficient  complexity  to 
properly  weigh  the  relative  merits  of  the  guidance 
variations    discussed  [Ref.    3]. 
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Figure  3, 1   Roland  Hissile  System  Operational  Schematic 
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17.     ON-OFF    (BANG^BANG)     CONTROL 

As  discussed  before,  LOS  guidance  maintains  a  missile 
position  on  the  LOS.  Usually  missile  position  has  a  cross 
range  error  (CRE)  and  we  want  to  reduce  this  error  to  zero 
in  the  minimum  time.  This  kind  of  problem  can  be  solved  by 
using  the   on-off  control.    The    basic   concept   of  this    is    that; 

Given  a  system  for  which  the  drive  is  limited  (has  a 
maximum  or  saturation  value) ,  the  fastest  response  is 
obtained  if  maximum  forward  drive  is  applied  at  t  =  0, 
and  is  reversed  at  'a  proper  instant  t  =  t1  so  that 
deceleration  under  maximum  reverse  drive  reduces  the 
velocity  to  zero  at  precisely  tha  command  value  of  the 
output.   The   drive   is  then   set   to    zero. 

Tha      ideal        relay      permits        only      two        conditions;        full 
acceleration  and   full  deceleration  [Bef.    5]. 

From  the  Bang-Bang  control  law,  we  can  derive  the 
switching  function  which  makes  the  error  go  to  zero  by  using 
the    proper    switching    time.       From    Newton's    second    law; 

F   =    m   CRE 

..            d           .                   F 
CRE    ■ (CRE)        = =      0 

dt  m 

CRE    =  CRE    dt         -        II    t    U1 

But    at  t   ■   0   ,    ere    ■   0  and    k1    =   0.    Therefore 

d 
CRE    = (CRE)      =   0    t  (4.  1) 

at 

CRE    =        CRE    dt      =   II    t2     +    k2  (4.2) 
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From    the   equation   4.1 

t   =      CRE  /    0 

t2    =     (CRE    /    U)2  (4.  3) 

Substitute      equation   4.3   into   equation   4.2 

CRE    =    (0   /    2)  (CRE   /    U)2      +    k2 

=     (CRE)2    /    (2    0)     ♦    k2  (4.4) 

where  lc2  is  integration  constant.  But  if  we  apply  a  full 
deceleration  at    the    halfway    point,    the   equation    4.4    becomes 

F    =     (CRE    |CRE|)     /(2    0)      ♦    CRE  (4.5) 

and   is   called  the  ERROR   FUNCTION.    U    will   be 
0    =     ±G 

or 

0    =    -(G)     SIGN(F)  (4.6) 

Equations  4.5  and  4.6  represent  the  SWITCHING  FUNCTION  which 
makes  the  error  go  to  zero  in  the  minimum  time.  The 
switching  function  and  the  block  diagram  of  the  on-off 
controller  are  depicted  on  the  Figures  4.1  and  4.2.  And  we 
can  obtain  the  cross  range  error,  CRE,  by  doubly  integrating 
U  with   the    initial    condition  of   CRE.    We    have 

CRE    =         0   dt    ♦    CRE(O) 

CRE    =         CRE    dt    +    CRE(O)  (4.7) 

The      simulation    results     of      these      equations   are      giver,     on 
Figures   4.3   through    4.7   and    the  computer   program    is   attached 
(see    Appendix   B)  . 
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Application  of 
Negative  Control 


Figure    4,1        Parabolic  Switching   Function 
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Figure  4.3        CRE   versus   Time 
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Figure  4.4        CBE   versus   Time 
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Figure   4.7        0  versus  Time 
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V-  BASIC  COMMAND  TO  LJNE^OF;  SIGHT  SIMULATION 

A.   SCENABIO 

The  engagement  was  designed  with  the  ground  tracker  and 
missile   launching  unit  located   at   tie    origin. 

The  target  was  flown  accross  the  first  quadrant  from  a 
position  4000  meters  on  the  x-axis  and  1000  meters  on  the 
y-axis  (4000,1000).  The  velocity  vector  of  the  target  was 
parallel  to  the  x-axis  and  magnitude  was  250  meters  per 
second. 

Since  most  missiles  need  a  few  seconds  of  boost,  the 
missile  is  not  controlled  during  this  time.  We  assumed  that 
the  missile  was  controlled  after  one  second  from  the  firing 
time  and  controlled  by  PROS  RAMMED  GUIDANCE  up  to  this  time. 
After  the  time  of  missile  "capture",  the  missile  was 
controlled  by  the  on-cffr  TVC  method  with  the  LOS  guidance 
law.    The  simplified    flow  chart   is   shown    in   the   Figure   5.1. 

In   order  to    simplify   the   problem,    we   assumed   that: 

1)  the  velocity  vector  of  missile  ,  Vm,  was  parallel  to 
the  LOS  between  the  target  and  origin  and  the  magni- 
tude of    Vm    was  constant,    503    meters   per   second; 

2)  the  LATAX  was  applied  to  the  missile  at  right  angles 
to  the  LOS.  This  was  a  reasonable  assumption  for 
this  kind  of  missile.  So  the  angle  <p  ♦  <<-  in  the 
Figure    2.2    is  almost   same   to    angle     <&  ; 

3)  the  measurement  noise  was  zero  so  we  omitted  the 
filter ,    F(s)  ; 

2 

4)  the  magnitude  of  LATAX  was  150  meters /second  which 
was   about   15    Gs.  ' 

The    geometry  depicted  in   Figure   5.2   summarizes  the   geometric 
situation. 
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Figure  5.1   Simplified  Plow  Chart  of  Basic  LOS  Command 
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Tracker 


Target 


Figure   5.2        Geoaetry   of   Basic  LOS   Guidance 

For   mathmatical      convenience    of    simulation,      we      need  to 
define   the   sign    of   the  CRE    and   the   LATAX   as   follow; 

♦1CRE|    :    When    the  missile   position   is    upper-side   of   LOS 
-|CRE|    :    when    the   missile    position   is   lower-side   of   LOS 
♦    I  U\      :    when    the  LATAX   is   upward   direction 
-    |  Q |      :    when    the  LATAX  is  dowaward    direction 
This    sign      was    based   on      the   positive        6^ which      is    defined 
when      6^  is   greater    than    dl  . 

B.       PBOGBAMHED    GUIDANCE    PHASE 


Since  the  major  emphasis  of  this  paper  was  on-off 
control,  we  assumed  that  the  missile  flew  along  the  LOS 
during  the  programmed  guidance  phass.  But,  in  a  practical 
situation,  there  is  some  cross  range  error  which  is  occured 
by      disturbances      such     as        wind,         propulsion      system      and 
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autopilot      time    delay,etc.  Hence      we    made      initialization 

errors,    and  the    on-off  control   started   with   these   errors. 

C.       ON-OFF,    THRUST    VECTOR,     MISSILE   CONTROL 

The  detail  of  the  on-off  control  was  discussed  before, 
hence  we  applied  this  to  the  LOS  guidance  scheme.  The  block 
diagram  of  this  system  is  depicted  in  the  Figure  5.3 
[Ref.    3]. 

In  order  to  determine  the  CRE,  the  tracker  estimates  the 
missile's  range  (Rm)  ,  by  the  elapsad  time  of  flight  and  the 
missile's  velocity  profile.  The  program  of  this  simulation 
is   attached  in    Appendix   C. 
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Figure  5.3    Block  Diagram  of  the  Basic  LOS  Coamand 


D.   SIMULATION  RESULTS 

Figure  5.4  shows  the  missile  and  target  geometry  in  X-Y 

plane.    The  missile  intercepted  the  target   at  the   point 

A(2605,10Q0)  with  the  almost  zero  miss  distance. 
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Figure  5.5(a)  shows  the  distance  between  target  and 
missile  versus  time.  The  distance  decreased  linearly  and 
neared   zero   at    the    time   at    5.58   seconds. 
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Pigure    5.4        The   Basic   LOS    Command 

Figure  5.5(b)  shows  the  CRE  varsus  time.  The  initial 
cross  range  arror  at  the  missile  capture  time  one  second 
after  firing  was  about  50  meters.  Since  tha  CRE  at  the 
"missile  capture  time"  was  positive,  the  CRE  increased 
initially.  So  the  maximum  CRE  was  about  58.2  meters  at  the 
time  1.330  seconds.  Then  it  decreased  to  almost  zero  meter 
at      2.55   seconds.  To   gat      a    fastar      response,       we      should 

increase  the      magnitude    of    the      LAT&X.       We    should      note   here 
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that  the  CRE  does  not  maintain  zero  value  because  we  did  not 
consider  the  target  motion  terms  in  this  phase  of  the  simu- 
lation [Bef-  2].  So  the  missile  had  some  small  cross  range 
error  and  the  BANG-BANG  controller  had  tried  to  reduce  this 
error  in   a   chatter-mode. 

Figure  5.6(a)  shows  the  CRE  versus  time.  Figure  5.6(b) 
shows  the  CRE  versus  CRE.  As  we  expected  this  curve  followed 
the    SWITCHING      FUNCTION   as      shown   in      figure    4.4.  Figures 

5.7(a)  and  5.7(b)  show  the  F  versus  time  and  the  0  versus 
time . 

This  program  was  tested  using  maneuvering  targets  and 
the  results  ware  almost  same  except  the  impact  position.  The 
results  cf  this  simulation  were  shown  on  the  Figures  5.8 
through  5.11  and  the  program  is  attached  in  Appendix  D.  The 
comparison   cf   these    simulations   is   summarized   in   Table   I. 


30 


2:2 

—  IX 


00 -a  00 'o        00 -a- 


U 


a 

at 

a 
m 
o 
u 

o 

-3 

u 

•H 
W 
(0 

CO 


en 

> 

0 

UJZ 

c_)cr 

zx 

err 

1—0 

tnc_) 

•— 1 

acj 

W 

•— < 

u 

>—  <r> 

"^ 

OCX 

S 

t— OD 

H 

1 

C/> 

_J 

H 

en 

a 

2: 

oo-oeri       aa'odn 


•oris       oo-ods 


00 "OS?         00*002 

,0t" 


a 

u 

r" 

fO 

4J 

en 

•w 

a 

CO 

e 

4J 

■H 

a> 

c-1 

cr 

- 

(0 

W 

H 

3 

1 

i/i 

cy 

M 

»H 

CU 

•  ^ 

> 

in 

in 

-C 


in 

M 

3 


31 


UJQ 
C_)d 

•  2: 

tno 

>(_) 

o 
o— 1 
ujcn 

IE  CX 


.S3 

u 


OO'O        oo-oz- 


00  -on-         00-09-         00-O»-       OO'OOI- 


T3 
O 
<Tj 

a 
s 
o 
u 

V) 

o 

-J 

u 

(D 

id 

ca 


L 


UJQ 
3TZ 

—  en 


en 

OCX 
UJCQ 

cc 


.S3 

u 


a 


00-08 


00 -on         M*£ 


cn 

a 

W 

M 

> 


.S3 


in 

M 

a 


32 


2:  ex 

1—3: 
o 

o">c_> 

>  — 
en 
<x 

XJCQ 


ooo       oo-st-     aoosi-     oo-k?-     ooooc- 


LU 

2:0 


(no 

>C_) 


.(_) 


S 


OO'K  boTS  OO'dK  00   i  WW  OO'il 


Cm 


p 

(0 

a 
■ 
o 
u 

in 

o 

- 

u 

•H 

0) 
<d 
QQ 

<I) 

J3 


in 

M 

a 


33 


>1 

u 
o 

■P 

u 

0) 
•n 
OS 

H 
H 

+J 
CD 

CP 

W 

i 

0) 

r-4 
•H 

cn 
1/1 

•H 

S 


4J 

0> 

or> 
W 
(0 
H 

en 

~ 

•H 
M 

> 

Q) 

O 

23 
0) 

J3 
£1 


•H 


oo-o7i       oo-ooi        oo 'do 
,0l» 


35        oVtP 


3U 


oo  ■  do         go -da         oo-d* oo-oa  oo-'o oo-o^- oo'oti-       oo-aa- oo-oi- oo-ooi'-° 


8 

r* 


u 


o 
cc 

a 


UJO 


c         oo'te         aa'di  do  s  ao-o         no's- 


U 


o 

M 
(0 

t^ 

a 

•H 

M 

a> 

> 
a 
<v 
a 
td 

13 
<D 


in 

- 

•H 

Bu 


35 


2 

H 
E-t 


UJ 

o 

(X 


— .O 


■  crz 
cduj 


—  a: 

<r 
ox 

QCI— 


ao-iJa         ao-as        oo'sn         oo-dn 


•3c         oo-te oo-di 


ao-9  oo'o         oo-e- 


en 


LJO 

en 

UJ 

•  > 

>UJ 

a: 
DX 
UJ 

EX 
t_)t— 


.22 


a  a  -so         oo'da         oo'dit         ooTa  «ro       oo-o^-       oo-oti-       oo'0»-       oo-oo-     oo'ooi- 


<J 


<P 

CP 
M 
(0 

cn 
fi 
•H 
M 
0) 
> 

a 

SC 


i      1 


0) 
U 

a 


36 


TIME 


a. qo  0.70 


Figure   5*11        0    versus    Time  for  a   Maneuvering   Target 


TABLE    I 
The  Basic    LOS   Command   Simulation   Result 


NO  N- MANEUVERING    TARGET 

.  time  (control)     =    1.0   sec 
CRE(O)  =    49.910(m) 

CRED(O)       =    49.832(m/sec) 

.time  (max. CRE)  1.33   sec 

=    53.  1 


CRE  (max) 


84(m) 


58    sec 


. time (intercept)  -  5, 
(Xm,Ym)    =     (2504.7,999.88) 
CRE  =    4.3344   E-5 

distance=    0.34894(m) 


MANEUVERING    TARGET 

1.0    sec 
49. 907  (m) 
49.891  (i/sec) 

1.33  (sec) 
58.201  (m) 

5.61    s  ec 
(2597.9,     1057.7) 
-6.5448    S-6 
0.47889(m) 
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71.     PSEUDO-LOS    COMMAND    SIMULATION 

The  guidance  scheme  cf  the  lead  angle  command  is  almost 
tha  same  as  that  of  the  basic  LOS  command.  Instead  of  the 
tracker-to-target  line-of- sight  this  guidance  scheme  uses 
the  tracker-to-estimated  impact  point  and  is  called 
"synthetic  line-of- sight"  (SLOS)  ,  or  "pseudo  line-of-sight". 
The  missile  is  controlled  to  fly  along  this  pseudo  line-of- 
sight.  The  block  diagram  of  this  system  is  easily  modified 
from    that   of  the  basic  LOS    and    is  shown    in    Figure   6. 1. 

The  estimated  impact  point  at  the  instantious  time  is 
calculated  by  using  the  "time  to  go"  (Tg)  and  the  "closing 
velocity"    (Vc>     between  the    target    and  the    missile. 
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Figure  6.1   Block  Diagra*  of  the  Pseudo  LOS  Command  System 
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The    "closing     velocity"   and    "time      to   go"    are      calculated  as 
fellow: 


.xJfc 


Vc  =       ((Vtx-Vmx)2  ♦    (Vty-Vmy)*) 

Tg  =    (distance    between   target   and   missile)    /   Vc 


f(Xt-Xm)2  +     (Yt-Ym)2)^    /   Vc 


The   missile      goes  to      the   impact:      point    directly.  The 

simulation  result  is  almost  same  as  in  the  basic  LOS  case. 
On  other  hand,  this  guidance  schema  is  poor  in  a  ECM  situ- 
ation. In  order  to  compare  the  results  we  used  the  same 
data  as  that  of  the  basic  LOS  command.  These  are  shown  in 
Figures  6.2  through  6.5  and  the  summarized  results  are  shown 
in   Table  II.    The  computer  program   is    attached   in    Appendix   E. 


1 

TABLE 

II 

The  Pseudo    LOS 

Co  mmand 

Simulation  Result 

.time  (control) 
CRE(O) 
CRED (0) 

1.0    sec 
49.912    (m) 
49.828     (m/sec) 

.time  (MAX.CRE) 
CRE  (max) 

1.33    sec 
58.  184     (m) 

.time  (intercept) 
(Xm,Ym) 

miss-distance 

5.58   sec 
(2604.7,    999.92) 
(2605.0,     1000.0) 
2.9572    E-5(m) 
0.35137     (m) 

_      i 
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Figure  6.5        0  versus  lime   for   the   Pseudo-LOS   Command 
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VII.     SIMULATIONS    WITH   TWO^LEVEL    RELAY    AND    SATURATION    CONTROL 

The  LOS  guidance  with  an  "ideal"  relay  has  beer. 
discussed.  In  this  chapter,  the  affect  of  the  different 
types      of  lateral     acceleration   demand      are   discussed.  In 

order  to  compare  the   results      with   the   previous   simulations, 
the    same   parameters    were  used. 

A.       TWO-LEVEL    RELAY 

The  large  magnitude  of  the  LATAX  makes  a  fast  response. 
But  in  "che  case  of  small  CREf  a  smaller  magnitude  of  LATAX 
is  needed.  This  idea  was  developed  in  a  "two-level1*  relay  as 
shown  on  the  Figures  7.1 (a;b) .  The  shaded  area  on  Figure 
7.2(b)  shows  the  region  of  a  lcwer  level  of  LATAX  in  the 
"CRE  verse  CRE"  phase  plain.  It  provided  the  minimum  over- 
correction. The  computer  programs  were  easily  modified  by 
adding   one   statement, 

IF  ((|CRE|  +  |3RSD|)  .LT.  M)  G  =  3/(N1/N2). 
We  used  the  values  150  m/sec2,  for  N1  and  15  m/sec2  for 
N2  and  1.0  for  M  in  the  simulations  of  the  basic  LOS  command 
and  the  pseudo  LOS  command.  The  results  were  almost  the  same 
as  the  previous,  except  in  the  figure  for  "U  versus  time". 
Table  III  summarized  these  simulation  results  and  Figures 
7.2,  7.3  and  7.4  show  the  "0  versus  time"  of  each  case.  The 
programs   were   attached  in   Appendix   F,    G    and  H. 
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Figure  7.1        Two-Level  Helay 


TABLE    III 
Two-Level  Relay  Control   Result 


.time  (control) 
CRE  (0) 
CRED(O) 

.time  (MAX.  CRE) 
CRE  (max) 

.  time  (intercept) 
Xm 
Ym 
Xt 
Yt 
CRE 
miss-distance 


NON-MVR 

MVR^TGT 

PSEUDO-LOS 

1.0 

1.0 

1.0 

49. 910 

49.907 

49.910 

49.  832 

49.891 

49.828 

1.33 

1.33 

1.33 

58. 184 

58.201 

58. 184 

)            5.58 

5.61 

5.58 

2604.7 

2597.9 

2604.7 

999.88 

1057.7 

999.92 

2605.0 

2597.5 

2605.0 

1000.0 

1057.5 

1000.0 

-2. 63E-6 

-3.63E-6 

5.03E-8 

0.34894 

0.47889 

0.35137 
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Figure   7.2        0*  versus   lime   for   the   Basic  LOS   Guidance 

B.       SATURATING    LINEAR  CONTROL 

In  the  previous  section  tha  two-level  relay  was 
discussed.  The  "sautrating  linear  control"  as  depicted  in 
the  Figure  7.5(a)  and  7.5(b)  was  also  studied.  The  shaded 
area  on  the  Figure  7.5(b)  shows  the  region  of  linear  control 
in  the  "CRE  versus  CRE"  phase  plans.  The  computer  programs 
were    easily  modified    by   adding    one   statement, 

IF     (ABS  (F)  .LE.  M)     U   =   -3    *    F   /   M. 

The  value  of  "H"  determines  the  linear  region  for  F. 
The  Figures  7.6  through  7.8  show  the  simulation  results  of 
the  basic  LOS  command  against  the  non-maneuvering  target 
case  for  "M"  equal  1,  5  and  10.  When  choosing  the  value  "M" 
equal  to  "one",  the  intercept  time  and  miss  distance  are 
almost  the  same  as  the  counterpart  of  the  ideal  relay  case. 
Hence   the  saturating   linear    control      can    be   used    in    practice 
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Figure  7.3        U    verse  Tiae  against  the   J1VR   Target    with 

Two-Level   Relay 
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Figure  7.4       a    verse  Time  of  the   Pseudo-LOS   with 

Two-Level  Relay 


U7 


instead  of  the  ideal  relay  by  choosing  a  proper  value  of 
"M".  The  summarized  results  are  in  the  Table  IV. 
Figures  7.9  through  7.11  show  the  results  of  the  maneuvering 
target  case  and  Figures  7.12  through  7.14  show  the  results 
of  the  pseudo-LOS  case.  These  programs  are  given  in 
Appendix   I,    J   and  K. 
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Figure   7.5       Linear   Switching   Relay 
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TABLE    IV 
Saturating  Linear  Control    Rasult    (H=1) 


UNIT 

•  time  (control) 
CRE  (0) 
CRED(O) 

.time  (MAX.     CRE) 
CRE  (max) 

. time (intercept) 
Xm 
Ym 
Xt 
Yt 
CRE 
miss-distance 


NON3MVR 

MVR^TGT 

PSEUDO-LOS 

1.0 

1.0 

1.0 

49. 910 

^9.907 

49. 912 

49.  832 

49.891 

49. 828 

1.33 

1.33 

1.33 

58. 184 

58.201 

58. 184 

5.58 

5.61 

5.58 

2604.7 

2598.0 

2604.7 

999.84 

1057.6 

999.88 

26  0  5.0 

2597.5 

2605.0 

1000.0 

1057.5 

1000.0 

-.0349 

-.0702 

-.0349 

0.3  507 

0.4841 

0.3492 
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fill.     CONCLUSION 

The  co m par i s io n  of  results  for  the  ideal  relay,  two- 
level  relay  and  saturating  linear  control  against  the  basic 
LOS  command  and  pseudo-LOS  command  against  a  non-maneuvering 
target  and  maneuvering  target  are  provided  in  the  Tables  V, 
VI  and  VII.  These  simulation  results  clearly  demonstrate 
that  "on-off"  control  of  a  missila  is  highly  desireable  and 
that  "saturating  linear  control"  of  a  missils  has  little 
adverse   effects    compared  to    an   "idaal   relay"    control. 


TABLE   7 
Comparison    of  the   Basic  LOS    Command  Results 


UNIT  IDEAL  TWO-LEVEL  SATURATION 
RELAY                       RELAY  CONTROL 

.timejccntrol)  1.0  1.0  1.0 

CRE(O)  49.910  49.910  49.910 

CRED(0)  49.832  49.832  49.832 

.time  (MAX.  CRE)  1.33  1.33  1.33 

CRE(max)  58.184  58.184  58.184 

.time  (intercept)  5.58  5.58  5.58 

Xm  2604.7  2604.7  2604.7 

Ym  999.88  999.88  999.88 

CRE  4.33E-5  4.33E-5  4.33E-5 

miss-distancs  0.34894  3.34894  0.34894 


59 


TABLE    VI 
Comparison   of  the  Maneuvering   Target   Results 


UNIT 


.time  (control) 
CRE(O) 
CRED(O) 

-time  (MAX.  CRE) 
CRE  (max) 

.time  (intercept) 
Xm 
Ym 
CRE 
miss-distance 


IDEAL 
RELAY 

1.0 

49. 907 
49.  891 

1.33 

58. 201 

5.6  1 
2597.5 
1057.5 
-6. 54E-6 
0.47889 


TWO- LEVEL 
RELAY 

1.0 

49.907 

49.891 

1.33 
58.201 

5.61 

2597.5 

1057.5 

3.63E-6 
3.47889 


SATURATION 
CONTROL 

1.0 

49.907 
49. 891 

1.33 

58. 201 

5.61 

2597.5 

1057.5 

-0.0702 

0.4841 


TABLE    711 

Comparison   of 

the  Pseudo- 

•LOS    Command 

Results 

UNIT 

IDEAL 
RELAY 

TWO- LEVEL 
RELAY 

SATURATION 
CONTROL 

.  time  (control) 
CRE  (0) 
CRED(0) 

1.0 

49. 912 
49.  828 

1.0 
49.910 

49.828 

1.0 

49.912 

49.828 

.time  (MAX.  CRE) 
CRE  (max) 

1.33 

58. 184 

1.33 

58.184 

1.33 

58. 184 

.time  (intercept) 
Xm 
Ym 
CRE 
miss-distanca 

5.58 
2604.7 

999.92 
2.96E-5 
0.35137 

5.58 
2604.7 

999.92 
5.03E-8 
0.35137 

5.58 
2604.7 

999.92 
2.96E-5 
0.35137 
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APPENDIX    A 
VARIABLES    LIST 


DIAGRAM 
VARIABLES 

COMPUTER 
VARIABLES 

CRE 

CRE 

CRE 

CRED 

F 

F 

G 

G 

0 

a 

6t 

SIGT 

6m 


SIGM 


dmt 

SIGMT 

vc 

VC 

tg 

TG 

t  (control) 

TCON 

SLOS 

SLOS 

Xm 

XM 

Ym 

YM 

Xt 

XT 

Yt 

YT 

Vm 

VM 

Vt 

VT 

Vmx 

VMX 

Vray 

VMY 

Vtx 

VTX 

Vty 

VTY 

NOON  DESCRIPTION 

cross -range- error 

rate  of   :ross-range-error   change 

error   function 

magnitude    of   lateral   acceleration 

missile's    lateral   acceleration 

angle    between   the   LOS    to   target 

and   X-axis 

angle    between   the    beam    to    missile 

and   X-axis 

angle  difference  between   m  and  t 

crossing  velocity 

time  to  go 

beginning    time    of   on-off    control 

synthetic    line-of-sight 

X-coordinate   of    missile   position 

Y-coordiaate   of    missile   position 

X-coordinate   of   target    position 

Y-coordinate   of    target    position 

velocity  of  the  missile 

velocity  of  the  target 

X-componsnt   of   missile's    velocity 

Y-componsnt  of   missile's    velocity 

X-component   of    target's   velocity 

Y-componsnt  of    target's   velocity 
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APPENDIX    B 
PROGRAM   OF    THE    SWICHING    FONCTION 


TITLE 

TITLE 

TITLE 

INTGER 

CONST 

INITIAL 


********** 

DERIVATIVE 
********** 

NOSORT 

* 


BANG-BANG    CONTROL 
SWITCHING       FONCTION 

*    YEUN,    J.Y.    * 
NPLd 
NPLOT=1 

CRE      =     1.0 
CRED    =    0. 
ACRED    =    0. 


SAMPLE 


TERMINAL 


G    =    1 

F    =   C 
0    =  - 
CRED 
CRE 
ACRED 
IF  (CR 

CALL 
CALL 
CALL 
CALL 
CALL 


.0 

RE  + 
G  *  SI 
=  INTG 
=  INTG 
=  ABS 
E.  LE. 


(CRED*ACRED)  /(2*G) 
GN(1.,F) 
BL(0.,U) 

BL(CRE,  CRED) 

(CRED) 

0.0)     CALL    ENDJ3B 


DRWG  (1  ,  1,  TIME,  CRE) 
, 1, CRE, CRED 
,1,TIME,F) 
,1, TIME, CRED) 


DRWG 
DRWG 
DRWG 
DRWG 


2 
3 
4 

5;i,time;u) 


CONT 
PRIN 
END 
STOP 


CALL    ENDRW(NPLOT) 
RL    FINTIM=2.  1,DELT=0.01  ,DELS  =  0.31 
T         0. 1 ,G,F,U, CRE, CRED 
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APPENDIX    C 
PROGRAM    OP   THE    BASIC    COMMAND    TO    LOS 


TITLE  BASI3    COMMAND   TO    L.O.S 

TITLE  WITH    IDEAL   RELAY 

TITLE  *****    YEUN,    J.Y.     ***** 

INTEG  RKSFX 

INTGER         NPLOT,KILL 

CONST  NPLOT=1 ,TCON=1.0 

CONST  VM=50  0. ,VT= 250, PI= 3. 141593,  KILL=0 

INITIAL 

XTO   =    U000. 

YTO   =    1 000. 

GAMT    =     PI 

F         =0. 
* 
***************************  ********************************* 

DERIVATIVE 

***************************  ********************************* 

NOSORT 
************************* 

***    TARGET    EARAMETERS    *** 
************************* 

VTX   =    VT*COS  (GAMT) 

VTY   '    VT*SINJGAMTf 
XT    =    VTX*TIME    ♦    XTO 
YT    =    VTY*TIME    +    YTO 
SIGT    =    ATAN2(YT,XT) 
IF     (IIME.GE.TCON)      GO    TO    50 
************************* 

***    MISSLE    PARAMETERS    *** 
************************* 

*  PROGRAMMED   GUIDANCE    ************************************** 
* 

SIGM    =    SIGT+0.  1 

VXM   =    VM    *    COS  (SIGM) 

VYM    =    VM    *     SIN(SI3M) 

XM      =    INTGRL(0.  ,VXM) 

YM      =    INTGRL(O..VYM) 

RM      =    SQRT(XM**2    ♦    YM**2) 

SIGMT    =    SIGM   -    SI3T 

CRE   =    RM    *    SIN  (SIGMT) 

CRED    =    DERIV(0.  ,CRE) 

GO    TO    200 
***************************  ** 

*  ON-OFF    GUIDANCE    (BANG-BANG    CONTROL)     ********************* 
***************************** 

* 

50  CONTINUE 

* 

G    =    150. 

ACRED    =    ABS(CRED) 
♦     fCRED*AC 


F    =   CRE    ♦     (CRED*ACRED)/(2*3) 
0    =   -G    *   SIGN(1.tFi 
CRED    =    INTGRLfO.O.U) 
CRE   =    INTGRL(CRE,CRED) 

RM      =    VM*TIME 
A3      =    CRE/RM 
SIGMT    =    ARSIN(A3) 
SIGM    =    SIGT    +    SI3MT 
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XM      =     RM    *     COS  (SIGM) 

YM       -    RM    *     SIN  (SIGM 
***************  ******************************************** 

200  CONTINUE 

* 
***************************  ******************************** 

***    MISSION    RESULT    ***      KILL    =    0     ;       TGT       HISSED 
**********************       KILL    =    1     ;       TGT       DESTROYED 

XDIST    =    XT-XM 

YDIST    =    YT-YM 

DIST       =    SQRT(XDIST**2    ♦    YDIST**2) 

IF     (DIST    .LE.5)        KILL    =     1 

IF    (DIST    .GT.5)        KILL    =    0 

IF     (XM     .GT.      (XT+30)  )  CALL    ENDJOB 
* 

SORT 

***************************  ******************************** 

********    OUTPUT    AND    PLOT    CONTROL    CARD    ********************* 
***************************  ******************************** 

SAMPL  E 

CALL  DRWG  (1  ,1,XM,YM) 

CALL  DRWG  (1,2,XT,YT) 

CALL  DRWG (2, 1f TIME, DIST) 

CALL  DRWG  (3,1, TIME, CRE) 

CALL  DRWG  (4,1, TIME, CRED) 

CALL  DRWG  (5  ,1,CRE,CRED) 

CALL  DRWG  (6,  1,TIME,F) 

CALL  DRWG  (7,1,  TIME  ,  U) 
TERMINAL 

CALL  ENDRW(NPLOT) 
CONTRL    FINTIM=6. 0,DELT=0.0  0 1, DEL 5=0.  0  03 
PRINT         0. 005, XM, YM, XT, YT, CRE, CRED, DIST, KILL 
END 
STOP 
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AP PENDIX  2 
PROGRAM  OF  THE  MANEUVERING  TARGET 


TITL 
TITL 
TITL 
INTE 
INTG 
CONS 
CONS 
INIT 


* 
**** 

DERI 
**** 

NOSO 
**** 
*** 
**** 


**** 

*** 

**** 

*   PR 

* 


E  BASIC    COMMAND    TO    LOS 

E  WITH    MANEUVERING    TGT 

E  ***      YEUN,     J.Y.       *** 

G  RKSFX 

ER         NPLOr,KILL 

T  NPLOT=1  ,TCON=1.0 

T  VM=50  0. ,VT=250,PI=3. 1415  93,  KILL=0 

IAL 

XTO   ■    4000. 

YTO   =    1000. 

GAMT    =     PI 

F    =   0. 

***********************  ********************************* 

VATIVE 
******************************************************** 

RT 
********************* 

TARGET    PARAMETERS    *** 
********************* 

VTX   =    VT*COS  (GAMT) 
VTY   =    VT*SINJGAMT) 
XT    =    VTX*TIME    ♦    XTO 
YT    =    100*SIN(0.5*PI*TIME)  +  !fT0 
SIGT    =    ATAN2(YT,XT) 
IF     (riME.GE.TCON)     GO    TO    50 
********************* 

MISSLE    PARAMETERS    *** 
********************* 

OGRAMMED   GUIDANCE    ************************************** 

SIGM    =     SIGT+0.1 

VXM    =     7M    *    COS  (SIGM) 

VYM    =    VM    *     SIN(SIGM 

XM      =    INTGRL(0.  ,7XM 

YM      =    INTGRL(0.  .VYMi 

RM      =    SQRT(XM**2    ♦    YM**2) 

SIGMT    =    SIGM   -    SIGT 

CRE   =    RM    *    SIN  (SIGMT) 

CRED    =    DERIV(0.  ,CRE) 

GO    TO    200 


***********************  ** 

-OFF    GUIDANCE    (BANG- BANG    CONTROL) 


**** 

*   ON 

***************************  ** 

* 

5 

* 


********************* 


0 


CONTINUE 


G    =    150 

ACRED    = 
F    =   CRE 


ABS  (CRED) 

♦     (CRED*ACRED)/(2*G) 
U    =    -G     *   SIGNM..F) 
CRED    =     INTGRLfO.O,  U) 
CRE   =    INTGRL(CRE,CRED) 
ACRE    =     ABS  (CRE) 

RM     ■=    VM*TIME 
A3      =    CRE/RM 
SIGMT    =    ARSIN(A3) 
SIGM    =    SIGT    ♦    SIGMT 
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XM      =     RM    *    COS  (SIGM) 
YM      ■    RM   *     SIN  (SIGM) 

:*******  ****************  ******************************** 

iO  CONTINUE 


****: 

20i 

***************************  ******************************** 

***    MISSION    RESULT    ***       KILL    =    0     ;       TGT       MISSED 
**********************      KILL    =    1     ;       TGI       DESTROYED 


MISSION    RESULT    ***       KILL    = 

:******************         KILL     - 

XDIST    =    XT-XM 

YDIST    =    YT-YM 

DIST      =    SQRT(XDIST**2 

IF     (DIST    .LE.5)        KILL 

IF    (DIST    .GI.5)       KILL    = 


* 

SORT 
**** 
**** 
**** 
SAMP 


IF     (XM 


YDIST**2) 

1 

0 


(XT  +  30)  )  CALL    ENDJOB 


**** 
**** 
**** 
LE 


*******************  ******************************** 

OUTPUT    AND    PLOT    CONTROL    CARD    ********************* 
*******************  ******************************** 


TERMINAL 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


DRWG 
DRWG 
DRWG 
DRWG 
DRWG 
DRWG 
DRWG 
DRWG 


1  ,1,  XM,  YM) 
1  ,2,XT,YT) 
2,1, TIME, DIST) 
3r  1,TIME,CRE) 
4,1, TIME, CRED) 
5,1,CRE,CRED) 
(6,  1,TIME,F) 
(7,1, TIME  ,U 


CONT 
PRIN 

END 
STOP 


CALL    ENDRW(NPLOT) 
RL    FINTIM=5. 9,  DELT=0.00 1,DELS=0.  0  03 

T  0.0  05, XM, Y M, XT, YT,C RE, CRED, DIST, KILL 
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APPENDIX    5 
PROGRAM    OP   THE    COMMAND    TO    PSEUDO-LOS 


TITLE  PSEUDO     -    LOS   COMMAND 

TITLE  WITH    IDEAL    RELAY 

TITLE  ***      YEUN,     J.Y.    ** * 

INTEG  RKSFX 

INTGER         NPLOT,KILL 

CONST  NPLOT=1  ,TCON=1.0 

CONST  VM=50  0. , VT= 250, PI= 3. 1415  93,  KILL=0 

INITIAL 

XTO   =    4000. 

YTO   =    1000. 

TG    =   0. 

f      =   0. 

GAMT    =    PI 
* 
***************************  ********************************* 

DERIVATIVE 
************************************************************ 

NOSORT 
************************* 

***    TARGET    PARAMETERS    *** 
***********  ***  *********** 

VTX   =    VT*COS(GAMr) 

VTY   =    VT*SINJGAMT) 

XT    =    VTX*TIME    «■    XTO 

YT    =    VTY*TIME    +    YTO 

SIGT    =    ATAN2(YT.XT) 

XLOS    =     XT    *    VTX*rG 

YLOS    =    YT    ♦     VTY*TG 

SLOS    =     ATAN2  (YLOS,  XLOS) 

IF     (TIME. GE. ICON)     GO    TO    50 
************************* 

***    MISSLE    PARAMETERS   *** 
************************* 

*  PROGRAMMED   GUIDANCE    ************************************** 
* 

SIGM    =    SLOS    +    0.1 

VMX   =    VM    *    COS(SIGM) 

VMY  =    VM   *    SINJSIGM 

XM      =    INTGRL(0.  ,VMX) 

YM      =    INrGRL(0.  .VMY) 

RM      =    SQRT(XM**2    *     YM**2) 

SIGMS    =    SIGM   -    SLOS 

CRE   =    RM    *     SIN(SIGMS) 

CRED   =    DERIV  (0.  fCRE) 

GO    TD    2  00 
***************************** 

*  ON-OFF    GUIDANCE    (BANG-BANG    CONTROL)     ********************* 
***************************  ** 

* 

50  CONTINUE 

* 

G    =    150. 

ACRED    =    ABS  (CRED) 
F    =   CRE    +     (CRED*ACRED)/(2*G) 
U    =   -G    *    SIGNM..F) 
CRED    =     INTGRLfO.O,  U) 
CRE   =    INrGRL(CRE,CRED) 
* 

RM      ■    VM*TIME 
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**** 

20 

* 

**** 
*** 
**** 


* 

SORT 
**** 

**** 
**** 
SAMP 


A3      =  CRE/RM 

SIGMS  =    ARSIN(A3) 
SIGM    =    SLOS    +    SIGMS 

VMX   =  VM    *     COS  (SLOS' 

VMY   =  VM    *     SIN(SL3S 

XM      =  RM    *    COS(SIGM 

YM      =  RM   *     SIN(SIGM' 

***********»*******?*************************  ***  ******* 

0  CONTINUE 

***********************  ******************************** 

MISSION    RESULT    ***      KILL    =    0     ;       TGT       MISSED 
******************      KILL    =    1     ;       TGT       DESTROYED 

XDIST    =    XT-XM 

YDIST    =    YT-YM 

DIST       =    SQRT(XDIST**2    +    YDIST**2) 

VC    =    SQRT (  (VTX-VMX) **2+ ( VTY- VMY) **2) 

TG    =    DIST/VC 

IF     (DIST    .LE.5)        KILL    =     1 

IF    (DIST    .GT.5)        KILL    =    0 

IF     (XM    -GT.      (XT+30)  )  CALL    ENDJOB 

***************  **************************************** 

****    OUTPUT    AND    PLOT    CONTROL    CARD    ********************* 
***********************  ******************************** 

LE 


TERMINAL 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


DRWG 
DRWG 
DRWG 
DRWG 
DRWG 
DRWG 
DRWG 
DRWG 


M*KM) 


1  ,2,XT,YT 
2/1, TIME, DIST) 
3,1, TIME, CRE) 
4,  1,TIME,CRED) 
5  ,1,CRE,CRED) 
6,1,TIME,F) 
7,1,  TIME  ,U) 


CONT 
PRIN 
END 
STOP 


CALL    ENDRWfNPLOT) 
RL    FINTIM=5. 9, DELT=0.00 1    ,DELS=0.003 
T  0.  1,TG,XM,YM  , XT, YT , CRE ,CRED, DIST , KILL 
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APPENDIX    F 
PROGRAH   OF   THE    BASIC   COMMAND    TO    LDS    WITH    TWO-LEVEL    EELAI 


TITL 
TITL 
TITL 
INTE 
INTG 
CONS 
CONS 
INIT 


* 
**** 

DERI 
**** 

NOSO 
**** 
*** 
**** 


E 

E 

E 

G 

ER 

T 

T 

IAL 


**** 

VATI 
**** 

RT 
**** 

TARG 
**** 


BASIC 

WITH 

RKSFX 
NPLOT 
NPLOT 
VM=50 

XTO   - 
YTO   = 
GAMT 
F      = 


COMMAND   TO    L.O.S 
TWO    LEVEL    RELAY 
YEUN,    J.Y.    ***** 


,KILL 

=  1  ,rcoN=i.o 

0. ,VT=250,PI=3.1415  93,KILL=0 

u  ooo. 


a  ooo 

1000. 

=  pi 
o. 


**** 

*** 

**** 

*    PR 


**** 

MISS 
**** 

OGRA 


*******************  ********************************* 

VE 
**************************************************** 

***************** 

ET    PARAMETERS   *** 
***************** 

VTX   =    VT*COS(GAMT) 
VTY   =    VT*SIN?GAMT) 
XT    =    VTX*TIME    +    XTO 
YT    =    VTY*TIME   ♦    YTO 
SIGT    =     ATAN2(YT/XD 
IF     (TIME. GE. ICON)     GO    TO    50 
********  ********* 

LE    PARAMETERS   *** 
***************** 

MMED   GUIDANCE    *****  **********************  *********** 


**** 

*   ON 
**** 

* 

5 

* 


SIGM 

VXM   = 

VYM   = 

XM      = 

YM      = 

RM      = 

SIGMT 

CRE   = 

CRED 

GO    TO 
********** 


=  SIG 
VM  * 
VM  * 
INTG 
INTG 
SQRT 
=  SI 
RM    * 

=  DER 
2  00 

***** 

HCE    £ 


-OFF    GOIDA 
***********************  ** 


T  +  0.1 

COS(SIGM) 

SIN(SIGM 
RL(0.,VXM 
RL  (0.  -VYM 
(XM**2    ♦     YM**2) 
GM    -    SIGT 

SIN  (SIGMT) 
IV  (0.  ,CRE) 

********  ** 

BANG-BANG    CONTROL) 


********************* 


CONTINUE 


G    =    150. 

ACRE    =    ABS(CRE) 

ACRED    =    ABS  (CRED) 

IF  ((ACRE+ACRED)     .LT.     1.)     G    =    15 

F    '   CRE    ♦     (CRED*ACRED)/(2*G) 

0    =   -G    *    SIGNM..F) 

CRED   =    INTGRL(0.0,O) 

CRE   =    INTGRL(CRE,CRED) 

RM      =    VM*TIME 
A3      =    CRE/RM 
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SIGMT    =    ARSIN(A3) 
SIGH    =    SIGT    ♦    SIGMT 
XM      =    RM    *    COS  (SISM) 
YM      =     RM    *    SIN(SIGM) 
*********** *******  *********  ****************  ******  ********** 

200  CONTINUE 

* 
***************************  ******************************** 

***    MISSION    RESULT    ***      KILL    =    0    ;       TGT       MISSED 
**********************      KILL    =    1     *       TGT       DESTROYED 

XDIST    =    XT-XM 

YDISr    =    YT-YM 

DIST      =    SQRT(XDIST**2    +    YDIST**2) 

IF    (DIST    .LE.5)        KILL    =    1 

IF     (DIST    .GT.5J        KILL    =    0 

IF     (XM     .ST.     (XT+30))CALL    ENDJOB 
* 

SORT 
*********************************************************** 

********    OUTPUT     AND    £LOT    CONTROL    CARD    ********************* 
***************************  ******************************** 

SAMPLE 

CALL  DRWG  (1  ,1,XM,YM) 

CALL  DRWG  (1  ,2,XT,YT) 

CALL  DRWG(2, 1,TIME ,DIST) 

CALL  DRWG    3,  1,TIME,CRS) 

CALL  DRWG (4,1, TIME ,CRED) 

CALL  DRWG  (5,1, CHE,  CRED) 

CALL  DRWG/6  ,1,TIME  ,F) 

CALL  DRWG  (7,  1, TIME, U) 
TERMINAL 

CALL  ENDRWjNPLOT) 
CONTRL    FINTIM=5.  9,DELT=0.00 1,DELS=0.  003 
PRINT         0.0  05,XM,YM,XT,YT,CRE,CRED,DIST,  KILL 
END 
STOP 
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APPENDIX    6 
PROGRAM    OF    THE    MANEUVERING    TARGET    WITH    TWO-LEVEL    RELAY 


TITL 
TITL 
TITL 
INTE 
INTG 
CONS 
CONS 
INIT 


* 
**** 

DERI 
**** 

NOSO 
**** 
*** 
**** 


E  BASIC    COMMAND   TO    LOS 

E  WITH    MANEUVERING    T  GT 

E  ***      YEUN,     J.Y.       *** 

G  RKSFX 

ER         NPLOT,KILL 

T  NPLOT=1,TCON=1.0 

T  VM=50  0. ,VT=250,PI=3.1415  93, KILL=0 

IAL 

XTO   =    4000. 

YTO   =    1000. 

GAMT    =    PI 

F      =    0. 

***************** ************************ *************** 

VATIVE 
******************************************************** 

RT 
********************* 

TARGET    PARAMETERS   *** 
********************* 

VTX   =    VT*COS(GAMT) 

VTY   =    VT*SINJGAMT) 

XT    =    VTX*TIME    +    XTO 

YT    =    100*SIN(0.5*PI*TIM2) +YT0 


**** 

*** 

**** 

*    PR 
* 


SIGT    =     ArAN2(YT.XT) 

IF     (TIME.GE.TCON)     GO    TO    50 


************  ********* 

MISSLE    PARAMETERS    *** 
********************* 

OGRAMMED   GUIDANCE    ************************************** 

SIGM    =    SIGT+0.1 

VXM   =    VM    *    COS  (SIGM) 

VYM   =    VM    *     SIN  (SIGM 

XM      -    INTGRL(0.  fVXM] 

YM      =    INTSRL  JO.  .VIM'i 

RM      =    SQRT(XM**2    «•     YM**2) 

SIGHT    =    SIGM    -    SIGT 

CRE   =    RM    *     SIN(SIGMT) 

CRED   =     DERIV(0.  ,CRE) 

GO    TO    2  00 
***********************  ** 

-OFF    GUIDANCE    (BANG-BANG    CONTROL) 


**** 

*   ON 
***************************  ** 

* 

50  CONTINUE 

* 

G    ■    150. 
ACRE    =    ABS(CRE) 
ACRED    =    ABS  (CRED) 
IF((ACRE*ACRED)  .LT.     1.)     G    =     15 
F    =   CRE    ♦     (CRED*ACRED)/ (2*3) 
U    ■   -G    *    SIGN(1..F). 
CRED    =     INTGRL(0.0,  U) 
CRE   =    INrGRL(CRErCRED) 
ACRE    =    ABS (CRE) 
* 

RM      =    VM*TIME 
A3       =    CRE/RM 


********************* 
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SIGHT    =    ARSIN(A3) 
SIGH   =     SIGT    +    SIGMT 
XM      =    RH    *    COS(SIGM) 
YM      =     RM    *     SIN  (SISM) 

***************   ******************************************** 

200  CONTINUE 

* 

***************************  ******************************** 

***    MISSION    RESULT    ***      KILL    =    0     ;       TGT       MISSED 
**********************      KILL    =    1     ;       TGT       DESTROYED 

XDIST    =    XT-XM 

YDIST    =    YT-YM 

DIST      =    SQRT(XDIST**2    ♦    YDIST**2) 

IF     (DIST    .  LE.5)        KILL    =    1 

IF     (DIST    .GT.5J        KILL    =    0 

IF     (XM     .ST.      (XT+30) ) CALL    ENDJOB 
♦  . 

SORT 
*********************************************************** 

********    OUTPUT     AND    PLOT    CONTROL    CARD    ********************* 
*********************************************************** 

SAMP  LE 

CALL  DRWG  (1  ,1,XM,YM) 

CALL  DRWG  (1  ,2,XT,YT) 

CALL  DRWG (2,1, TIME , DIST) 

CALL  DRWG  (3,1,  TIME,  CRE) 

CALL  DRWG(4, 1, TIME  , CRED) 

CALL  DRWG  (5,  1, CRE,  CRED) 

CALL  DRWG(6,1,TIME  ,F) 

CALL  DRWG  (7,  1, TIME, U) 
TERMINAL 

CALL  ENDRW  (NPLOT) 
CONTRL    FINTIM=5.  9,DELT=0.00 1,DELS=0.  003 
PRINT         0.005, XM,YM, XT, YT,  CRE, CRED,  DIST,  KILL 
END 
STOP 
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APPENDIX    H 
PROGRAM    OF   THE    COMMAND    TO    PSEUDO-LOS    WITH    TWO-LEVEL    RELAY 


TITLE  PSEUDO    -    LCS    COMMAND 

TITLE  WITH    TWO-LEVEL    RELAY 

TITLE  *****    YEUN,    J.Y.     ***** 

INTEG  RKSFX 

INTGER         NPLOr,KILL 

CONST  NPLOT=1  ,TCON=1.0 

CONST  VM=50  0. ,7T  = 250, PI  = 3. 1 41 5  93,  KILL=0 

INITIAL 

XTO   =    4  000. 

YTO   =    1000. 

TG    =    0. 

GAMT    =    PI 

F      =    0. 
* 
***************************  **********************  *********** 

DERIVATIVE 

***************************  ********************************* 

NOSORT 
************************* 

***    TARGET    PARAMETERS    *** 
************************* 

VTX   =    VT*COS(GAMT) 

VTY   =    VT*SINJGAMIJ 

XT    =    VTX*TIME   +    XTO 

IT    =    VTY*TIME   ♦    YTO 

SIGT    =    ArAN2(YT.XT) 

XLOS    =     XT    ♦    VTX*TG 

YLOS    =     YT    ♦    VTY*TG 

SLOS    =    ATAN2(YLOS,  XLOS) 

IF    (HME.GE.TCON)     GO    TO    50 
************************* 

***    MISSLE    FARAMEIERS    *** 
************************* 

*  PROGRAMMED   GUIDANCE    ************************************** 
* 

SIGM   =    SLOS    ♦    0.1 

VMX   =    VM    *    COS  (SIGM) 

VMY   =    VM    *    SIN  (SIGM 

XM      =    INTGRL(0.  ,VMX 

YM      =    INTGRLJO..VMY 

RM      =    SQRT(XM**2    ♦    YM**2) 

SIGMS    =    SIGM    -    SLOS 

CRE   =    RM    *     SIN(SIGMS) 

CRED    =    DERIV(0.  ,CRE) 

GO    TO    200 
***************************  ** 

*  ON-OFF   GUIDANCE    (BANG-BANG    CONTROL)     ********************* 
***************************  ** 

* 

50  CONTINUE 

* 

G    =    150. 

ACRE    =    ABS(CRE) 

ACRED    =    ABS(CRED) 

IF  ((ACRE*  ACRED)     .LT.     1.)     G    ■    15. 

F    =   CRE    ♦     (CRED*ACRED)/ (2*3) 


U    =    -G    *   SIGN(1..F) 
CRED    =    INTGRLIO.O,  U) 
CRE   ■    INTGRL(CRE,CRED 


) 
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RM      =  VM*TIME 

A3      =  CRE/RM 

SIGMS  =    ARSIN(A3) 

SIGH   =    SLOS    ♦    SIGMS 

VMX   =  VM    *    COS(SLOS] 

VMY  =  VM    *    SINfSLOS 

XM      =  RH   *    COS  (SIGM 

YM      =    RM    *    SIN(SIGM) 
***************************  ******************************** 

200  CONTINUE 

* 
*********************************************************** 

***    MISSION    RESULT    ***      KILL    =    0     ;       TGT       MISSED 
**********************      KILL    =    1     ;       TGT       DESTROYED 

XDIST    =    XT-XM 
YDISI    =    YT-YM 

DIST      =    SQRT(XDIST**2    +    YDIST**2) 
VC    =    SQRT ( (VTX-VMX) ** 2+ ( VTY- VMY) **2) 
TG    =    DIST/VC 

IF     (DIST    .LE.5)       KILL    =     1 
IF     (DIST    .GT.5)        KILL    =    0 
IF     (XM    .GT.      (XT  +  30))CALL    ENDJOB 
* 

SORT 

***************************  ******************************** 

********    OUTPUT    AND    PLOT    CONTROL    CARD    ********************* 
*********************************************************** 

SAMPLE 

CALL  DRWG  (1  ,  1,XM,YM) 

CALL  DRWG(1  ,2,XT,YT) 

CALL  DRWG (2, 1,TIME, DIST) 

CALL  DRWG  (3,1, TIME, CRE) 

CALL  DRWG(4,  1,TIME,CRED) 

CALL  DRWG(5  ,1,CRE,CRED) 

CALL  DRWG  (6,1,TIME,F) 

CALL  DRWG  (7,1,  TIME  ,  U) 
TERMINAL 

CALL  ENDRW(NPLOT) 
CONTRL    FINTIM=5. 7,DELT=0.00 1    ,DELS=0.003 
PRINT          0.  1,TG, XM, YM  , XT, YT, CRE ,C RED, DIST , KILL 
END 
STOP 
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APPENDIX    I 
PROGRAM    OF    THE    BASIC    COMMAND   TO    LOS    WITH    SATURATION    CONTROL 


TITLE  BASI3    COMMAND    TO    L.O.S 

TITLE  MISSLE      CONTROL 

TITLE  WITH    SATORATION    CONTROL 

INTEG  RKSFX 

INTGER         NPLOIr KILL- CUR 

CONST  NPLOT=1 ,TCO N-1 . 0,C UR= 1 

CONST  VM=5Q0. ,VT= 250, PI= 3. 1 41 5  93,  KILL=0 

PARAM  M    =    10. 

INITIAL 

XTO   =    4000. 

YTO   =    1000. 

GAMT    =    PI 

F      =    0. 
* 
***************************  ********************************* 

DERIVATIVE 

***************************  ********************************* 

NOSORT 
************************* 

***    TARGET    PARAMETERS    *** 
************************* 

VTX   =    VT*COS(GAMT) 
VTY   =    VT*SINJGAMT| 
XT    =    VTX*TIME    +    XTO 
YT    =    VTY*TIME   +    YTO 
SIGT    =     ATAN2(YT,XI) 
IF     (TIME. GE. ICON)     GO    TO    50 
****************  ********* 

***    MISSLE    PARAMETERS   *** 
************************* 

*    PROGRAMMED   3UIDANCE    ************************************** 
* 

SIGM    =     SIGT+0.1 

VXM   =    VM    *     COSfSIGM] 

VYM   =    VM    *     SIN  (SIGM 

XM      =    INTGRLfO.  ,VXM 


YM      =    INrGRL(0.  .VYM} 

RM      =    SQRT(XM**2    ♦     YM**2) 

SIGMr    ■    SIGM   -    SIGT 

CRE   =    RM    *     SIN  (SI3MT) 

CRED    ■     DERIV(0.  ,CRE) 

GO    TO    2  00 
*** * **************  ********* ** 

*   ON-OFF    GUIDANCE    (BANG-BANG    CONTROL)     ********************* 
**** **************** ******* ** 

* 

50  CONTINUE 

* 

G    =   150. 

ACRED    =    ABS  (CRED) 

F    =   CRE    ♦     (CRED*ACRED)/(2*G) 

U    =   -G    *    SIGN(1.,F) 


IF     (ABS(F)      .LT.    M)     0=    -G*F/M 
CRED    =    INTGRL(0.0,U) 
CRE   =    INTGRL(CRE,CRED) 


RM      =    VM*TIME 
A3      =    CRE/RM 
SIGMT    ■    ARSIN(A3) 
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SIGM    =    SIGT    +    SIGMT 
XM       =    RM    *     COS  (SIGM) 
YM      =     RM    *    SINfSIGMJ 

*************#-****##******* ******************************** 

200  CONTINUE 

***************************  ******************************** 

***    MISSION    RESULT    ***      KILL    =    0     ;       TGT       MISSED 
**********************      KILL    =    1     ;       TGT       DESTROYED 

XDIST    =    XT-XM 

YDISr    =    YT-YM 

DIST      =    SQRT(XDIST**2    ♦    YDIST**2) 

IF    (DIST    .LE.5)       KILL    =    1 

IF     (DIST    .GT.5)        KILL    =    0 

SORT 

***************************  ******************************** 

********    OUTPUT    AND    PLOT    CONTROL    CARD    ********************* 
***************************  ******************************** 

SAMPLE 

CALL    DRWG  (1  , CUR, XM.YM) 

CALL    DRWG(2  , CUR, TIME,  DIST) 

CALL    DRWG  (3  , CUR, TIME,  CRE) 

CALL    DRWG (4  ,CUR,CRE,CRED) 

CALL    DRWG  (5  , CUR, TIME,  F) 

CALL    DRWG  (6  , CUR, TIME,  U 
TERMINAL 

IF     (3UR    .EQ.    3)     CALL    ENDRW(NPLOT) 

CUR   =    CUR    ♦     1 
CONTRL    FINTIM=5.  65  ,DELT=0.  0  0  1 ,  DELS=3  .  00  3 
PRINT         0.005,XM,YM, XT, YT, CRE, CRED, DIST, KILL 
END 

PARAM  M   =    5. 

END 

PARAM  M    =    1. 

END 
STOP 
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APPENDIX    J 
PROGRAfl    OF    THE    MANEUVERING    TARGET    IITH    SATURATION    CONTROL 

TITLE  BASIC    COMMAND    TO    LOS 

TITLE  (MANEUVERING    TGT) 

TITLE  WITH    SATURATION    CONTROL 

INTEG  RKSFX 

INTGER         NPLOr, KILL-CUR 

CONST  NPLOT=1 ,1  CO N=1 . 0, CUR= 1 

CONST  VM=50  0. ,VT=250, PI  a  3 . 1 41 5  93 , KILL=0 

PARAM  M    =    10. 

INITIAL 

XTO   =    4000. 

YT0  =    1000. 

GAMT    =     PI 

F      =    0. 
* 

***********  A*************************************  *********** 

DERIVATIVE 

***************************  ********************************* 

NOSORT 
************************* 

***    TARGET    PARAMETERS    *** 
************************* 

VTX   =    VT*COS  (GAMT) 
VTY   =    VT*SINjGAMTJ 
XT    =    VTX*TIME    ♦    XTO 
YT    =    100*SIN(0.5*PI*TIME)  +?T0 
SIGT    =    ATAN2JYT,XT) 
IF     (TIME.  GE. ICON)      GO    TO    50 
************************* 

***    MISSLE    PARAMETERS    *** 
**************  *********** 

*  PROGRAMMED  GUIDANCE    ************************************** 
* 

SIGM    =     SIGT+0.1 
VXM   =    VM    *    COS  (SIGM] 
VYM   =    VM    *     SIN  (SIGM 
XM      =    INrGRL(0.,VXM 
YM      =    INTGRL(0.  .VYM) 
RM      =    SQRT(XM**2    ♦    YM**2) 
SIGMT    =    SIGM   -    SIGT 
CRE  =    RM    *    SIN  (SIGMT) 
CRED    =    DERIV(0.  ,CRE) 
GO   TO    2  00 
***************************  ** 

*  ON-OFF    GUIDANCE    (BANG-BANG    CONTR3L)     ********************* 
***************************  ** 

* 

50  CONTINUE 


G    =    150 
ACRED    =    AB 

U    = 


=    ABS  (CRED) 
RE    ♦     (CRED*ACRED)/(2*G) 
G    *   SIGN(1.,F) 


IF     (ABS  (F)  . LT.M)     U=-G*F/M 

=     INTGRLfO.O,  U) 
CRE   =    INTGRL(CRE,CRED) 


CRED    =     INTGRLfO.O 
CRE   =    INTGRL(CR~ 
ACRE    =    ABS  (CRE) 


RM      =    VM*TIME 
A3      =    CRE/RM 
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SIGHT    =    ARSIN(A3) 
SIGM   =    SIGT    «■    SIGMT 
XM       ■    RM    *    COS  (SIGM) 
YM      =    RM    *    SIN(SIGM) 

*******#*#*******#**#*##****£#*******#*#****#*$**  ********** 

200  CONTINUE 

* 

***************************  ******************************** 

***    MISSION    RESULT    ***      KILL    =    0     ;       TGT       MISSED 
**********************      KILL    =    1     ;       TGT       DESTROYED 

XDIST    =    XT-XM 

YDIST    ■    YT-YM 

DIST      =    SQRT(XDIST**2    ♦    YDISI**2) 

IF    (DIST    .  LE.5)       KILL    =    1 

IF     (DIST    .GT.5)       KILL    =    0 
* 

SORT 

***************************  ******************************** 

********    OUTPUT    AND    PLOT    CONTROL    CARD    ********************* 
*********************************************************** 

SAMP  LE 

CALL    DRWG  (1  ,CUR,XM,YM) 

CALL    DRWG  (2  ,CUR,TI  ME, DIST) 

CALL    DRWG  (3, CUR,  TIME,  CRE) 

CALL    DRWG (4  ,CUR,CRE,CRED) 

CALL    DRWG  (5, CUR, TIME,  F) 

CALL    DRWG  (6, CUR, TIME,  U) 
TERMINAL 

IF     (CUR    .EQ.3)     CALL    ENDRW(NPLOT) 

CUR    =     CUR    +     1 
CONTRL    FINTIM=5. 65  ,  DELT=0.  0  01 ,  DELS=0  . 00  3 
PRINT         0.005, XM,YM, XT, YT, CRE, CRED, DIST, KILL 
END 

PARAM       M    =    5. 
END 

PARAM       M    =    1. 
END 
STOP 
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PROGRAM    OF    THE    COMMAND    TO    PSEUDO-LOS    WITH    SATURATION    CONTROL 


TITLE 
TITLE 
TITLE 
INTEG 
INTGE 
CONST 
CONST 
PARAM 
INITI 


R 


AL 


* 
***** 

DERIV 
***** 

NOSOR 
***** 
***  T 
***** 


***** 
***  m 
***** 

*    PRO 

* 


PSEODO    -    LOS   COMMAND 

WITH      SATURATION    CONTROL 

*****    YEUN,    J.Y.    ***** 

RKSFX 

NPLOT, KILL, CUR 

NPLOT=1  ,TCON=1.0,CUR=1 

VM=50  0. ,VT=250,PI=3.1415  93,KILL=0 

M    =   10. 

XTO   =    4000. 
YTO   =    1000. 
TG   =   0. 
GAMT    =    PI 
F      =    0. 

******************************************************* 

ATIVE 

**********************  ********************************* 

T 
******************** 

ARGET    PARAMETERS   *** 
******************** 

VTX   =    VT*COS(GAMT) 
VTY  =    VT*SINJGAMT) 
XT   =    VTX*TIME   +    XTO 
YT    =    VTY*TIME    ♦    YTO 
SIGT    =     AIAN2(YT,XT) 
XLOS    =     XT    ♦    VTX*T3 
YLOS    =     YT    ♦    VTY*TG 
SLOS    =    ATA N2j YLOS,  XLOS) 
IF    (TIME.  GE.TCON)     GO    TO    50 
******************** 

ISSLE    PARAMETERS    *** 
******************** 

GRAMMED  GUIDANCE    ************************************** 


SIGM   =    SLOS    ♦    0.1 
VMX   =    VM    *    COS  (SIGM) 
VMY   =    VM    *    SIN(SI3  M 
XM      =    INTGRL 
YM      =    INTG 
RM      =    SQRT 
SIGMS    =    SIGM   -    SLOS 
CRE   =    RM    *    SIN  (SIGMS) 
CRED    =    DERIV(0.  ,CR  E) 


[     SIN (513  M) 

;rl(0.  ,vmx 

IRL(0.  .VM  Y' 
'(XM**2    ♦    Y 


M**2) 


GO    TO    2  00 
**********************  ** 

OFF   GUIDANCE    (BANG-BANG    CONTR3L 


***** 

*   ON- 
***************************** 

* 

50  CONTINUE 

* 


********************* 


G    =    150. 

ACRED    =    ABS(CRED) 

F    =   CRE    ♦     (CRED*ACRED)/(2*3) 

U    =   -G    *   SIGN(1.,F) 

IF     (ABS(F)      -LT.M)      U   =    -G    *    F    /M 

CRED    =     INTGRLfO.O,  U) 

CRE   =    INTGRL(CRE,CRED) 
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* 

RM      =  VM*TIME 

A3      =  CRE/RM 

SIGHS  =    ARSIN(A3) 

SIGM    =     SLOS    +    SI3MS 

VMX   =  VM    *    COS  (SLOS] 

VMY   =  VM    *     SIN(SLOS 

XM      =  RM    *    COS  (SIGM 

YM      =    RM    *     SIN  (SIGM 
**** * ************* *********  **********************  ********** 

200  CONTINUE 

* 
*********************************************************** 

***    MISSION    RESULT    ***      KILL    =    0     ;       TGT       MISSED 
**********************      KILL    =    1     ;       TGT       DESTROYED 

XDIST    ■    XT-XM 

YDIST    =    YT-YM 

DIST      =    SQRT(XDIST**2    ♦    YDIST**2) 

VC    =    SQRTMVTX-VMX)  **2+ (VTY-VMY)  **2) 

TG    =    DIST/VC 

IF    (DIST    .LE.5)       KILL    =    1 

IF     (DIST    .GT.5)       KILL    =    0 
* 

SORT 

****************  ***  ********  **********************  ********** 

********    OUTPUT    AND     PLOT    CONTROL    CARD    ********************* 
***************************  ******************************** 

SAMPL  E 

CALL    DRWG  (1  ,CUR,XM,YM) 

CALL    DRWG(2  , CUR  ,TI  MEr  DIST) 

CALL    DRWG  (3  , CUR, TIME,  CRE) 

CALL    DRWG (4  ,CUR ,CR E ,CRSD) 

CALL    DRWG(5, CUR, TIME,  F) 

CALL    DRWG  (6, CUR, TIME,  U) 
TERMINAL 

IF     (CUR    .EQ.3)     CALL    ENDRW(NPLOT) 

CUR   =     CUR    ♦     1 
CONTRL    FINTIM=5. 65 ,DELT=0.  0 01     ,DELS=0.003 
PRINT         0.005,TG,XM, YM, XT, YT, CRE, CRED, DIST,  KILL 
END 

PARAM  M    ~    5. 

END 

PARAM         H   ■    1. 
END 
STOP 
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